Special stability advantages of position-Verlet over velocity-Verlet in multiple-time step integration
نویسندگان
چکیده
We present an analysis for a simple two-component harmonic oscillator that compares the use of position-Verlet to velocity-Verlet for multiple-time step integration. The numerical stability analysis based on the impulse-Verlet splitting shows that position-Verlet has enhanced stability, in terms of the largest allowable time step, for cases where an ample separation of time scales exists. Numerical investigations confirm the advantages of the position-Verlet scheme when used for the fastest time scales of the system. Applications to a biomolecule, a solvated protein, for both Newtonian and Langevin dynamics echo these trends over large outer time-step regimes. © 2001 American Institute of Physics. @DOI: 10.1063/1.1389855#
منابع مشابه
Optimized particle-mesh EwaldÕmultiple-time step integration for molecular dynamics simulations
We develop an efficient multiple time step ~MTS! force splitting scheme for biological applications in the AMBER program in the context of the particle-mesh Ewald ~PME! algorithm. Our method applies a symmetric Trotter factorization of the Liouville operator based on the position-Verlet scheme to Newtonian and Langevin dynamics. Following a brief review of the MTS and PME algorithms, we discuss...
متن کاملEfficient multiple-time-step integrators with distance-based force splitting for particle-mesh-Ewald molecular dynamics simulations
We develop an efficient multiple-time-step force splitting scheme for particle-mesh-Ewald molecular dynamics simulations. Our method exploits smooth switch functions effectively to regulate direct and reciprocal space terms for the electrostatic interactions. The reciprocal term with the near field contributions removed is assigned to the slow class; the van der Waals and regulated particle-mes...
متن کاملVerlet-i/r-respa Is Limited by Nonlinear Instability
Abstract. This paper uncovers additional stability limitations of multiple time stepping (MTS) integrators for molecular dynamics (MD) that attempt to bridge time scales. In particular, it is shown that when constant-energy (NVE) simulations of Newton’s equations of motion are attempted using the MTS integrator Verlet-I/r-RESPA, there are nonlinear instabilities when the longest step size is a ...
متن کاملTime-reversible quantum molecular dynamics
Time-reversible ab initio molecular dynamics based on a lossless multi-channel decomposition for the integration of the electronic degrees of freedom [Phys. Rev. Lett. 97, 123001 (2006)] is explored. We present a lossless time-reversible density matrix molecular dynamics scheme. This approach often allows for stable Hartree-Fock simulations using only one single self-consistent field cycle per ...
متن کاملPractical 03: Position-Based Dynamics (e.g., Verlet System)
This practical focuses on position-based dynamics using the Verlet integration scheme. The student needs to implement a real-time interactive soft-body effect using Verlet mechanics (e.g., cloth with collision detection, forces, and movement).
متن کامل